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In high temperature QCD, the perturbation theory is plagued with infrared di- 
vergences which reflect long-range non-perturbative phenomena. I argue that it is 
possible to study such phenomena within a classical thermal field theory which can 
be put on a three-dimensional lattice. The classical theory is an effective theory 
for the soft, non-perturbative modes, as obtained after integrating out the hard 
modes in perturbation theory. It is well suited for numerical studies of the non- 
perturbative real-time dynamics, which cannot be studied within the standard, 
imaginary-time formulations of lattice QCD. 



1 Introduction 

There is currently a considerable interest in the physics of ultrarelativistic plas- 
mas, by which we mean plasmas in, or near, thermal equilibrium at very high 
temperature, much larger than any mass scale: T ^> to. This interest, which 
has triggered many theoretical advances, is mainly motivated by two impor- 
tant applications: The first is to the deconfined phase of QCD, the celebrated 
quark-gluon plasma (QGP), which is expected to be found in the heavy ion 
collisions at RHIC and LHC. The second is to the high-temperature, symmet- 
ric, phase of the electroweak theory, which is relevant for baryogenesis in the 
early Universe. 

If the temperature is high enough, the physics of non-Abelian gauge theo- 
ries is expected to be simple. This common wisdom, which we shall see below 
to be a little too optimistic, is motivated by asymptotic freedom: as the tem- 
perature increases, the coupling "constant" becomes small, g(T) -C 1, and 
perturbation theory becomes applicable. That is, the high-T plasma can be 
treated as a gas of weakly interacting quarks and gluons []. 

By using perturbation theory, significant progress has been indeed achieved 
in understanding the long-wavelength excitations of the plasma and the related 
screening phenomena (see Refs. an for a summary and more references, and 
Sec. 2 below for.-a brief account). Furthermore, the free energy in hot QCD has 
been computed 13 up to the order g 5 , which is the highest accuracy permitted 

6 For definitness, I shall mostly use a QCD-inspired terminology. Note, however, that most of 
the present considerations do also apply to the electroweak plasma, provided the temperature 
is high enough (T > T c ). 
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by perturbation theory (the correction 0(g 6 ) turns out to be non-perturbative; 
see Sec. 3 below). 

At the same time, lattice simulations have been used to compute thermo- 
dynamic quantities □ (like the free energy and the^ntropy) and to study static 
aspects of the deconfininga and the electroweak □ phase transitions. Besides 
providing "exact" numerical results, the lattice calculations allow us to explore 
non-perturbative physics like the strong coupling regime in QCD at interme- 
diate temperatures (T ^ T c ~ 200 MeV). As the temperature increases, the 
comparaison with perturbation theory becomes possible: it has been found, 
for instance, that the lattice estimates for the energy density approach rather 
quickly the Stefan-Boltzmann limit above T c , thus comforting the picture of 
hot QCD as a weakly interacting gas []. 

It thus may come as a surprise that, in some cases, lattice simulations 
cannot be avoided not even in the study of the high temperature limit, where 
the coupling constant is arbitrarily small. This is so since perturbation the- 
ory is plagued with infrared divergences which reflect large collective effects 
(cf. Sec. 3). Moreover, such divergences affect not only static characteristics, 
like the free energy, but also dynamical quantities, that is, real-time correla- 
tion functions, for which the standard lattice simulations — as formulated in 
imaginary time — are not applicable. An important example, to which Lshall 
return later, is the rate of baryon number violation at high temperatureO. To 
compute such quantities, we need new, non-perturbative, methods which allow 
for numerical studies of the real-time dynamics. It is my purpose in this talk 
to present such a method which is based upon a semiclassical approximation 
(cf. Sec. 4). 

2 Collective behaviour and screening 

I consider a purely Yang-Mills plasma ^| in thermal equilibrium at a tempera- 
ture T which is high enough for the coupling constant to be small: g -C 1. In 
the absence of interactions, this would be simply the "black body radiation" , 
that is, a collection of free, massless gluons with typical energies of the order T 
and Bose-Einstein occupation numbers: Nq{E) = l/{c^ E — 1), with (3 = 1/T. 
However, the gauge interactions — although weak — do significantly change 
this picture. They give rise to an hierarchy of scales: 

T > gT > g 2 T > g A T ■ ■ • , (1) 

c Rather surprisingly, it turns out that the perturbative expansion for the free energy E is 
poorly convergent up to temperatures as high as 3 GeV. Still, as shown in Refs. Q, the 
convergence can be greatly improved by using Pade approximants. 
d kt high-T, quarks are not important for the infrared physics to be discussed here. 
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with the various scales corresponding to different physical phenomena. 

Thus, the typical excitations of the plasma are "hard" gluons, with mo- 
menta k ~ T. Such gluons can develop a collective behaviour over a typical 
space-time scale ~ 1/ gT, which is large as compared to the mean interparticle 
distance ~ l/T. This results in long-wavelength (A ~ ^s> 1/T) oscilla- 

tions of the averge colour density which are most economically described in 
terms of kinetic equations El. In these equations, the hard (k ~ T) gluons are 
represented by their average colour density SN a (k, x) to which couple the soft 
(i.e., long- wavelength) colour fields A£(x). The relevant equations read: 

d 3 fc 

(v-D x ) ab SN b (k,x) =- 9 vE.(a)-^, (2) 

where = + igA£T a is the covariant derivative, E l a = is the chromo- 
electric field, C A = N for SU(JV), and w M = (1, v) with v = k/fc denoting the 
velocity of the hard particles (k = |k|, so that |v| = 1). 

The first equation above is the Yang-Mills equation for the soft fields A%. 
It involves, in its r.h.s., the colour current induced by the collective motion of 
the hard particles: 



(D v F Vfl ) a (x) = 2gC A I v»6N a (k,x), 



d 3 fc 



jS(x) = 2gC A I —3 v» SN a (k, x). (3) 



In turn, the soft colour wave A£(x) acts as a driving force for the collective 
behaviour, and this is described by the second Eq. (Q) (whictumay be seen as 
a non- Abelian generalization of the familiar Vlasov equation Era) . By solving 
this equation, we can express the current fa m terms of the gauge fields A% 
and thus obtain an effective Yang-Mills equation which involves the soft fields 
alone: 

D v F^ = ml (4) 
D J 4tt v ■ D y ' 

Here, the angular integral J df2 runs over the orientations of the unit vector 
v, and mo is the Debye mass : 







Eq. (Q) describes the propagation of long- wavelength colour waves in the high- 
T plasma. The induced current in its r.h.s. is the result of the wave scattering 
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off the hard thermal particles. In general, this current is non-local, and also 
non-linear in the gauge fields (note the covariant derivative in the denomina- 
tor). Still, for time-independent fields, it reduces to a very simple expression: 
j^ (x) = (5 M °m|jA°(x), which defines a screening "mass" for the Coulomb field 
A° a . Indeed, for such static fields, the [i = component of Eq. (||) simplifies 
to: 

D •E + to^A (x) = 0, (6) 

which is the non-Abelian generalization of the Poisson-Debye equation: 

(-A + m£)i4o(x)=0, (7) 

and implies the screening of any electrostatic fluctuation A$ (x) over distances 
r ~ 1/mo- In other terms, the static (fco ~~ * 0) limit of the Coulomb propagator 
— as following from Eq. (^) — reads: 

A)o(fc = 0,k) = ^t-^t, (8) 

so that the Debye mass acts as an infrared cutoff ~ gT in the electric sector. 

The situation in the magnetic sector is more complex: in the static limit the 
vector current j a vanishes, as alluded to before, so that the time-independent 
magnetic fields are not screened. The analogue of Eq. (||) reads then: 

Dij (k = 0,k) = Sij ~J lk] , (9) 

which is the same as at tree level. For time- dependent magnetic fields, however, 
screening does.-occur, through the mechanism of Landau damping. This too 
can be studiedM on Eq. (||), with the result that, for small but non- vanishing 
frequency fc , the magnetic piece of the gluon propagator reads: 



- » - < 10 > 



Note the purely imaginary character of the self-energy in the denominator: 
this is a dissipative effect, describing the absorbtion of the magnetic field by 
the plasma constituents. Eq. (|o|) also shows that, for large enough frequencies 
feo ~ k, the magnetic fields are screened as efficiently as the electric ones. 

The above kinetic-theory picture of the screening (which is actually equiv- 
alent to the one-loop picture in the "hard thermal loop" approximation Era) 
turns out to be further complicated by non-perturbative phenomena. There 
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are indeed theoretical arguments 13, which are also supported by lattice sim- 
ulations and which suggest that screening should occur also for the static 
magnetic fields, but only at the softer scale g 2 T (see Sec. 4 below). In prac- 
tical calculations, this is often parametrized by introducing an infrared cutoff 
/i ~ g 2 T in the magnetic sector ("magnetic mass"). But it is fair to say that 
the corresponding physical mechanism is not yet fully understood. This is so 
since, as we shall see shortly, g 2 T is precisely the scale where perturbation 
theory breaks down. 

3 The breakdown of the perturbation theory 

The screening phenomena greatly improve the infrared behaviour of the per- 
turbation theory (PT), thus allowing for many interesting calculations! 2 ]. Still, 
when going to higher orders in PT, one is often confronted with severe infrared 
(IR) divergences, which signal a breakdown of the perturbation theory at the 
scale g 2 T. There are essentially two reasons for that: 

The first is the singularity in the magnetostatic propagator as k = |k| — > 
(cf. Eq. (P)). Although this might be screened, as alluded to before, at the scale 
g 2 T, nevertheless such a screening will not be enough to restore perturbation 
theory (see below). 

The second reason is the Bose-Einstein amplification of the soft modes 
with momenta k <C T: such modes have large thermal occupation numbers, 



and therefore give large "radiative" corrections in higher orders. Specifically, 
by adding a new loop to a preexistent Feynman graph, we generate a correction 
of relative order g 2 No(k), where k is the momentum carried by the added gluon 
propagator. If k ~ g 2 T, then N (k) ~ T/k ~ 1/g 2 and g 2 N (k) ~ 1 : that is, 
by adding more soft (k ~ g 2 T) loops, we remain at the same order in g, and 
the loop expansion breaks down. n 

To see a specific example of this difficulty, let's follow LindeEl and consider 
the higher order corrections to the free energy: a typical n-loop (with n > 4) 
diagram has the "ladder" topology in Fig. [I] and gives a contribution ~ 
g 6 T 4 (<7 2 T//z)™ (this is a simple power counting estimate). In this equation, 
fi is an infrared cutoff which has been introduced by hand to give a meaning 
to an otherwise IR-divergent loop integral. As we know by now, such a cutoff 
is indeed generated via the screening effects. In the electric sector, we have 
Debye screening and therefore /i ~ vtld ~ gT (cf. Eq. @); with ^ ~ gT, 
F( n > ~ g™ +2 T 4 , and higher loops contribute to higher orders in g, as it should 



N (k) 



I 




(11) 
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Figure 1: n-loop ladder diagram contributing to the free energy. 



for PT to make sense. In the magnetic sector, on the other hand, we have 
at most /x ~ g 2 T ', in which case all the diagrams with four or more loops 
contribute to the same order in g (namely, to the order g 6 ). We thus face a 
breakdown of PT in the magnetic sector which, unlike the electric sector, is 
not protected by Debye screening. 

For static quantities like the free energy, the non-perturbative corrections 
can be estimated, at least in principle, via lattice QCD. But there are also 
time- dependent correlation functions which appear to be non-perturbative and 
for which the standard lattice calculations (as formulated in imaginary time) 
are not applicable. Let me give you some examples in this sense: 

In relation to baryogenesis, one is interested in the rate for anomalous 
baryoiunumber violation in the high-T, symmetric, phase of the electroweak 
theory!! This is a genuinely non-perturbative phenomenon where the variation 
AB of the baryon number is tied up — via the chiral anomaly — to the 
transitions between topologically inequivalent vacua: 

AB(t) <xJ Q dx J d 3 xF- v Pr, (12) 

with F£ v = (l/2)e A " ypA Fp A . At high temperature, such transitions will neces- 
sarily involve very soft (k ~ g 2 T) magnetic fields, the only ones to be respond 
sible for non-perturbative phenomena. Still, it has been recently recognized □ 
that the topological transitions are also sensitive to the hard (k ~ T) plasma 
modes, via the Landau damping alluded to before (cf. Eq. (M)). 

But non-perturbative effects are also met in the study of simpler corre- 
lation functions, like the gluon, or quark, propagator: e.g., when studying 
the spectrum of the elementary excitations (the plasma "quasiparticles"), we 
consider the large time behaviour of the corresponding 2-point functions. On 
general grounds, one expects an exponential decay of the quasiparticles, as 
coming from their scattering off the plasma constituents. For instance, if 
S(t, p) = < ^(t, p)i/)(t, — p) > denotes the quark propagataor, then one ex- 
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Figure 2: Fermion-fermion elastic scattering in the Born approximation. The bubble on the 
gluon line indicates the resummation of the screening effects. 



pects S(t — > oo, p) oc e 1 " B ^'e~ 7 ' ! '' f 1 where E(p) is the mass-shell energy and 
7(p) is the damping rate, i.e., the_total interaction rate of the quark in the 
plasma. Still, explicit calculationslI9 of j(p) show a logarithmic infrared diver- 
gence already in leading order, as coming from collisions with the exchange of 
soft magnetic gluons (cf. Fig. |^). Similar divergences occur for gluons, and 
also for electrons in a hot QED plasma. Thus, the lifetime of the quasiparticles 
turns out not to be computable in perturbation theory. 

4 Classical effective theory for real-time processes 

I now present a method which allows us, at least in principle, to compute 
the non-perturbative, real-time n correlation functions alluded to before. 

The basic idea is not new 0: since the non-perturbative phenomena are 
associated with soft (k ~ g 2 T) magnetic fields which have large thermal oc- 
cupation numbers (cf. Eq. (p"l|)), some semi-classical approximation should be 
applicable. To see this in a simple way, let's put back Planck's constant h in 
Eq. ( |il"| ) and compute the average energy per mode in thermal equilibrium: 

hk 

£ {k) = ~ T as hk <C T. (13) 

As h — > 0, we recover the classical equipartition theorem, as expected. But 
the above example shows that the relevant inequality is hk -C T, so that the 
classical limit (h — > at fixed k and T) is actually equivalent to the soft 
momentum limit (fc — > at fixed h and T). This observation is useful since 
we know how to perform real-time lattice simulations for a classical thermal 
field theory: All we have to do is to solve the classical equations of motion for 
given initial conditions, and then average over the classical phase space with 
the Boltzmann weight exp(— (3H). Since the initial conditions (say 0(x) and 
0(x) for a scalar theory) depend only on the spatial coordinate x, the phase 
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space integration is actually a three-dimensional functional integral, which can 
be implemented on a lattice in the standard way (and actually with even less 
numerical effort than in four dimensions !). The only question is, what is the 
correct classical theory ? 

It has been originally assumedQ that, in order to compute the hot baryon 
number violation (cf. Eq. (E2|)), it should be enough to consider the classical 
Yang-Mills theory at finite—temperature. This hypothesis, which led to the 
lattice calculations in Ref. E3, is highly non-trivial: it assumes that the topo- 
logical transitions are totally insensitive to the hard (k ~ T) plasma modes, 
which are not properly described by the classical theory (since the approxima- 
tion in Eq. ( |l3|) fails at momenta k ~ T). And indeed, the classical theory is 
well-known to run into ultraviolet (UV) problems, like the famous "ultravio- 
let catastrophe" of Rayleigh and Jeans: The classical estimate for the energy 
density of the black body radiation, namely (A is an ad-hoc UV cutoff): 

W < i4 > 

is obviously wrong (since UV divergent), in contrast to the quantum result: 

which is finite since the large momenta k T arc exponentially suppressed 
by the Bose-Einstcin distribution function. 

Now, as already argued in Sec. 3, the topological transitions in hot QCD 
are driven by soft (k ~ g 2 T) field configurations, and thisJs_jthe main jus- 
tification for using the classical Yang-Mills theory in Refs. DO Still, it has 
been also argued in Sec. 2 that the dynamics of the soft modes is strongly 
modified by the hard particles, which generate screening. Thus, even though 
a soft process, the baryon number violation might still be sensitive to the hard 
particles, via the screening effects. And actually we have both theoretical □ 
and numerically evidence that this is indeed the case. 

Thus, in order to properly compute the non-perturbative correlation func- 
tions of interest, we need to correct the classical Yang-Mills theory by including 
the screening effects. That is, the relevant classical theory should be an effec- 
tive theory which applies only to the soft modes, but where the hard modes 
have been integrated over to generate screening. From Sec. 2, we have a 
candidate for such a theory: Eq. (^) includes indeed the screening effects, via 
the colour current in its r.h.s. Still, there are a few "technical" complications 
associated with this equation: First, Eq. (|J) is non-local (and also dissipative: 
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recall the imaginary part in the denominator of Eq. (|10[)), so it is not a priori 
clear how to construct the thermal phase-space and the Hamiltonian. Sec- 
ond, in order to avoid overcounting, we need a precise, and gauge-invariant, 
separation between hard and soft degrees of freedom. 

The discussion in Sec. 2 suggests a solutionEJ to the first problem above: 
rather than working with the non-local equation ^, we can conveniently replace 
it with the coupled system of local equations in Eq. (||). There is a price to be 
payed for that: in addition to the gauge fields A£(x), the local description in 
Eq. (g) also involves the average colour density SN a (k, x) , which can be seen 
as an "auxiliary field" . Still, when working with a local theory, we are in a 
better position to look for a Hamiltonian formulation, as I discuss now. 

The first step is to recognize, on the second Eq. (||), that the v and k- 
dependence can be factorized in SN a (k, x) by writing: 

5N a (k, x) = -gW a (x, v) (dN Q /dk). (16) 

The new functions W a (x,v) satisfy the equation: 

{v- D x ) ab W b (x,v) = vE„(a;), (17) 

which is independent of k since the hard particles move at the speed of light: 
|v| = 1. By using Eq. (|l6|), the induced current can be written as: 

f a tx) = m 2 D J^v»W a {x,v), (18) 

where the radial integration (i.e., the integration over k = |kj) has been ex- 
plicitly worked out, and the Debye mass mo is defined iHiEq. (]5|). 

The Hamiltonian formulation of the effective theoryEjO involves the aux- 
iliary fields W a (x, v) together with the soft gauge fields A£(x). In the temporal 
gauge Aq = 0, the independent degrees of freedom are Ef, A" and W a , and 
the corresponding equations of motion read ^: 

K = -d a A«, 

-d E? + <,j,il) J lh,:'' =mlJ^ Vl W a (x, v), 

(do + v • D) afc Wb — v • E a , (19) 
together with Gauss' law which in this gauge must be imposed as a constraint: 

G°(x) = (D • E) a + m 2 D J ^ W a (x, v) = 0. (20) 

e Notc that Eqs. (|is|) are not in canonical form: this is already obvious from the fact that 
we have an odd number of equations. 
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Eqs. ([L9|) are conservative; the cnrresponding, conserved energy functional 
(which also acts as a HamiltonianEil) has the following simple form El: 

H = \J d 3 z |E a -E a + B a -B a + m 2 D J ^ W a (x, v) W a (x, v)| , (21) 

which is manifestly gauge invariant. 

We are now in position to write down the classical partition function and 
compute (generally time-dependent) thermal expectation values. As discussed 
at the end of Sec. 3, we are interested in correlation functions of the magnetic 
fields A l a . These can be obtained from the following generating functional: 

Z cl [J9] = J V£?VA a i VW a 8{Q a ) exp j-/3H + J d 4 x J? 

where A l a (x) is the solution to Eqs. fll9| ) with the initial conditions {£?, Af , W a } 
(that is, Ef(to,x) = (x), etc., with arbitrary t ), and Q a and Ji are expressed 
in terms of the initial-fields (cf. Eqs. ( pp| ) and (pi])) 



(^(x) ,(22) 



It can be verifiedtil that the phase-space measure VE^VAfDW 1 in Eq. (22) 
is invariant under the time evolution described by eqs. (|l9|), so that Z c i[J] is 
independent of the (arbitrary) initial time to, as it should. (This point is not 
trivial because of the non-canonical structure of the equations of motion.) 

There is another essential — but technically quite involved — point that 
I am currently glossing over: this is the intermediate cutoff separating hard 
from soft degrees of freedom, and which should appear as an ultraviolet cutoff 
in Eq. ( |2^ ) (without such a cutoff, the effective theory would develop linear 
UV divergences to one loop order). As discussed in Ref. Eil, this cutoff can 
be indeed introduced in such a way to make the effective theory UV finite, 
but cancel — in the calculation of physical quantities — against appropriate 
"counter-terms" in the Hamiltonian. Moreover, the cutoff procedure proposed 
in Ref. Ed can be also implemented on a lattice; this is important since it allows 
one to take the continuum limit in the lattice calculations. 

For illustration, let me finally consider two simple, yet non-trivial, appli- 
cations of the effective theory. The first is the Jf = limit of Eq. ( p2| ) which 
yields, after some simple algebra, 



J VA%VA a t exp |-| J d 3 x (B?Bf + (XV4„) 



mj,A a A a ,(23) 



where the Aq components of the gauge fields have been reintroduced as La- 
grange multipliers to enforce Gauss' law. This is "almost" the thermal partition 
function for classical Yang-Mills theory: the only new feature is the screening 
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Figure 3: A generic self-energy correction which yields IR divergences in perturbation theory. 
All the diagrams of this ki nd a re resummed in the non-perturbative calculation of the fermion 
propagator based on Eq. (|22[), 



mass for the electric fields, which is the only trace of the screening effects in 
the static limit (cf. Eq. (||)). 

The expression in Eq_ig^) also coincides with the first order result of the 
"dimensional reduction" ali^L-a method which consists in integrating out the 
non-static Matsubara modes u to obtain an effective 3-dimensional theory for 
the static one. By putting this theory on a lattice, on£ has been able to perform 
accurate studies of the electroweak phase transitional and also to compute the 
non-perturbative contributions to the free-energy E3 (cf. Sec. 3) and to the 
non-Abelian Debye masslij. 

Note also that the magnetic sector of Eq. ( p3| ) (that is, the sector involving 
only the vector fields Af(x)) is formally the same as 3-dimensional Euclidean 
QCD with dimensionfull coupling constant 173 = gyT, This theory is ill be- 
haved in perturbation theory (the IR divergences in Fig. 1 may be seen, as 
an illustration of such a bad IR behaviour B), but it is generally believed £3 to 
generate a dynamical mass gap oc g| = g 2 T, which is the only mass scale in 
the problem. This is the celebrated "magnetic mass" alluded to in Sec. 3. 

As a second application, consider the large-time behaviour of the electron, 
or quark, propagator, which we have seen to be ill defined in perturbation 
theory (cf. Sec. 3). Fermions have not been yet included in the effective the- 
ory, but it is easy to compute the fermion propagator in a soft background 
field A 1 (in the eikonal approximations) , and then average over the thermal 
fluctuations of the background field as shown in Eq. (^). In QED, all these 
calculations can be done explicitly (the corresponding functional integral is 
Gaussian), with the striking result that^t very large times, the electron prop- 
agator shows a non- exponential decayo: S(t) oc exp{— aTt ln(mui)}, with 
a = e 2 /47r. In terms of diagrams, this calculation corresponds to a resumma- 
tion of all the "quenched" self-energy corrections illustrated in Fig. 3. 
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In QCD, the corresponding calculation can be performed only numerically, 
which requires a lattice implementation of the effective theory. More generally, 
such an implementation would allow for systematic studies of the real-time 
non-perturbative dynamics in hot gauge theories. 
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